Optimal Entangling Capacity of Dynamical Processes 
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We investigate the entangling capacity of dynamical operations when provided with local ancilla. A compari- 
son is made between the entangling capacity with and without the assistance of prior entanglement. An analytic 
solution is found for the log-negativity entangling capacity of two-qubit gates, which equals the entanglement 
of the Choi matrix isomorphic to the unitary operator. Surprisingly, the availability of prior entanglement does 
not affect this result; a property we call resource independence of the entangling capacity. We prove several use- 
ful upper-bounds on the entangling capacity that hold for general qudit dynamical operations, and for a whole 
family of entanglement monotones including log-negativity and log-robustness. The log-robustness entangling 
capacity is shown to be resource independent for general dynamics. We provide numerical results supporting a 
conjecture that the log-negativity entangling capacity is resource independence for all two-qudit unitaries. 
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Quantum information theory, and its various practical ap- 
plications, has lead to a mature theory of entanglement as a re- 
source. It is now common to think of entanglement as almost a 
fungible commodity. Two different quantum states with mani- 
festly very different non-locality properties may sometimes be 
converted into each other by local operations, they can also be 
diluted Jill and distilled 131 ■ The laws governing the capabili- 
ties of local operations to interchange between different entan- 
gled states are understood in terms of entanglement measures 
and monotones that quantify the amount of entanglement in a 
given quantum state pl] . 

Entanglement can be created between physical systems 
where there is a suitable non-local dynamical process. As 
we can quantify the entanglement of quantum states, we 
can quantify the efficacy of a dynamical operation at pro- 
ducing entanglement. Such quantities are called the entan- 
gling power, capacity or strength |]3l of a dynamical process. 
The average entanglement generated by unitaries, commonly 
called the entangling power, has been studied by Zanardi and 
coworkers ISJ-ull- Whereas the maximum entanglement that 
can be produced by unitaries is usually called the entangling 
capacity Isl- UOll . In the practical context of attempting to 
maximize entanglement production for utilization as part of 
a quantum information protocol, the most appropriate quan- 
tity is the entanglement capacity (herein the EC). 

Kraus and Cirac considered the EC for two-qubit unitaries 
acting on two-qubit product states f^. Leifer et at extended 
this analysis to find the entangling capacity when the initial 
two-qubit state has prior entanglement [^l, and found that this 
can boost the entangling capacity. In both papers it was ob- 
served that access to local ancilla can enable much higher en- 
tangling capacities for some unitaries, with the swap gate as 
the most striking example. Without local ancilla the swap gate 
has zero entangling capacity, but with local ancilla the swap 
gate can produce two Bell pairs. Despite numerous examples 
where local ancilla prove beneficial, no previous work has es- 
tablished an analytic solution for the entangling capacity with 
access to local ancilla. 

Furthermore, optimising the entanglement produced from a 
dynamical process, with access to local ancilla, is a problem 
that occurs in a broad range of architectures for quantum tech- 
nologies lUll - USll . Figure[T]shows the essential components of 
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FIG. 1: Many architectures for quantum communication and dis- 
tributed quantum computation have the above structure, possibly 
repeatedly over many locations. Physical systems A and B can- 
not freely interact (e.g. because of spatial separation). How- 
ever, they can interact by a fixed dynamical process described by a 
completely-positive trace-preserving map, A, which acts on Hilbert 
space Ha®Hb- Furthermore, physical systems A and B also contain 
quantum systems Ha' and Hy respectively. For example, both A and 
B may contain many identical qubits, such as in two ion traps con- 
taining many atoms I13il . We show finite ancillary Hilbert space, but 
for simplicity assume arbitrarily large local Hilbert spaces. Further- 
more, we assume unitaries and measurements within each location 
can be implemented with arbitrary precision. 



a physical system where this optimisation is applicable. For 
example, this structure is embodied by two separate eletro- 
magnetic traps each containing many ions, with a fixed optical 
process for producing entanglement between traps |13]]. 

We prove that two qubit unitaries have an EC, measured 
by the log-negativity 11141 USll . that has a simple closed form. 
Interestingly this EC equals the entanglement of the Choi ma- 
trix pu isomorphic to U via the Choi-Jamiolkowski isomor- 
phism jig,! 1711 . This proof will also hold for a restricted class 
of higher dimensional unitaries that we characterize. Cirac 
et at ilSll have previously proposed quantifying capabilities 
of entangling unitaries by considering the entanglement of 
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as measured by: (thick dashed) the von 



Neumann entropy with initially separable states; (thick solid) the 
von Neumann entropy assisted by prior entanglement; (thin) the log- 
negativity both with and without access to prior entanglement. 



the Choi matrix. However, Cirac et al did not show that the 
Choi matrix captures the maximum achievable entanglement 
for any continuous monotones. Our results give a concrete 
operational meaning to the entanglement of the Choi matrix. 
It is surprising that our result for two-qubit unitaries holds 
with or without the assistance of prior entanglement. We 
call this property resource independence of the log-negativity 
EC for two-qubits. In deriving our result for two-qubit uni- 
taries, we prove two theorems that are useful in their own 
right. These theorems provide upper-bounds on the EC in the 
more general setting of two-qudit completely-positive trace- 
preserving maps. The first upper-bound we prove also ap- 
plies to a whole family of EC monotones that includes both 
log-negativity EC and log-robustness EC. Our upper-bounds 
entail that log-robustness EC is always resource independent, 
and give us reason to suspect that log-negativity EC is also 
generally resource independent. We support this conjecture 
with a numerical study of the log-negativity EC for two-qutrit 
unitaries, where we found no violations of resource indepen- 
dence. 



I. NOTATION AND DEFINITIONS 

We aim to solve the following problem: for a given 
completely-positive trace-preserving (CPTP) map A acting on 
Hilbert spaces Ha ® Ht, what is the maximum entanglement 
that can be produced? This maximization problem depends 
on what resources are available. Throughout, we assume that 
local ancilla (typically labelled as Hilbert space Ha' (8) H^' ), 
local unitaries and classical communication are freely avail- 
able as is appropriate in many communication and distributed 
computation scenarios. The maximum achievable entangle- 
ment for a dynamical map. A, when the initial state is separa- 
ble, is given by: 



where S denotes the set of separable density matrices, and 
E{p) is an appropriate entanglement monotone for density 
matrices. However, when prior entanglement is available, the 
maximum increase of entanglement is: 



S^(A) = sup{^(A(a)) - E{a)\a G V}, 



(2) 



where V is the set of all physical density matrices, and we 
have denoted the availability of this resource by a superscript. 
This quantity cannot be smaller than the unassisted capacity, 
E{A), but in general may be larger. The von Neumann en- 
tropy 12511 . EyN, is an example where some unitary gates have 
been shown to exhibit an increase, i?^([/) > EyN{U), when 
assisted by entanglement lll9ll . We do not need to invoke ex- 
otic unitary gates to observe the phenomena, as it occurs even 
for simple two-qubit phase gates, as shown in Fig.|2] Alterna- 
tively, we may be interested in knowing how much entangle- 
ment could be produced given access to a restricted class of 
resources, such as bound entanglement 112011 : 



E^{A) = sup{E{A{a)) - E{a)\(T e B}, 



(3) 



where B is the set of all physical density matrices with only 
bound entanglement. Formally, all operators in B are positive 
semi-definite before and after partial transposition. 

Since the von Neumann EC violates resource indepen- 
dence, such that i?^ 7^ Ey]^, there is a temptation to infer 
that this phenomena would be observed for other EC mono- 
tones. Here we present results that counter this intuition, 
showing that many other entanglement monotones have an in- 
herent resource independence, which makes them cleaner, and 
simpler, candidates for EC monotones. 

We are principally interested in the logarithmic versions 
of negativity 11141 llSfl and robustness ll21ll and other similar 
decomposition-based monotones (defined in section|IIl). The 
negativity is defined as: 



En{p)^{\\p^\\~1)/2 

where F denotes a partial transpose over system B, and 
denotes the trace norm: 



= tr f \/IU 



(4) 



(5) 



The robustness of p was first defined Kill as the amount of 
any separable state a' that must be mixed with p to make the 
whole mixture separable: 



Es{p) = mi{tW 



ta' 



1 + t 



a, a' eS}, 



(6) 



where S is again the set of separable states. Despite the con- 
trast between the conventional definitions of E^ and Es, we 
will see that they are actually very closely related 11 1411 . 

Throughout, we are interested in logarithmically rescaled 
variants of both monotones as E^^ip) = log2(l + 2Ex{p))- 
Logarithmic variants of negativity and robustness are strongly 
subadditive under tensor products, such that: 



E{A) ^ sup{E{Aia))\a e S}, 



(1) 



E{pi®p2) <E{pi) + E{p2). 



(7) 



We desire this property as any monotone violating this condi- 
tion can never give rise to an EC that is resource independent 
(see AppendixlAj. Indeed, the concept of resource assisted EC 
can become meaningless when subadditivity is not respected. 
Conversely, if an EC is resource independent, then a rescaling 
of the underlying entanglement monotone that destroys sub- 
additivity must also remove resource independence. Under 
fairly weak assumptions such a rescaling can always be found 
(see AppendixlBJ). 

When a subadditive entanglement monotone always satu- 
rates the inequality, such that 



Ar='D^ , cT^>0 



Eln{pi ® /02) = Eln{pi) + Eln{P2), 



(8) 



the entanglement monotone is said to be strongly additive. 
Strong additivity gives entanglement monotones a clearer op- 
erational interpretation, and so is typically desirable. Any EC 
monotone that is both resource independent and strongly ad- 
ditive has the following elegant property: given a dynamical 
map A with entangling capacity, £'(A), we can use the dy- 
namical map n times to produce, at best, exactly nE{K) en- 
tanglement. The log-negativity is both strongly additive and 
provides an upperbound on the number of singlets that can be 
distilled 12211 . and so we will focus on this monotone. 

Note that, our notation for EC monotones always uses a cal- 
ligraphic superscript to denote the resources that are relevant 
to the maximization problem, and a subscript to describe the 
underlying entanglement monotone for quantum states. 



II. DECOMPOSITION-BASED MONOTONES 

Here we introduce the idea of decomposition-based mono- 
tones of entanglement, as first discovered by Werner and Vi- 
dal III4I1 . They proposed entanglement monotones of the form: 

Em{p) =ini{t\p= {l+t)a+-t(j--t>Q-(j^ C,M}. (9) 

This formula returns the minimum value of t over all real lin- 
ear decompositions of p into matrices ct^ that belong to a 
specific set M.. We will call operators in the set A^ zero- 
entanglement operators as for all cr e A^ we have Em (c) = 
0. However, remember that a zero-entanglement operator is 
not necessarily separable, as the monotone may fail to detect 
entanglement for some states. In order to ensure that the quan- 
tity is well-defined and monotonic under local operations, we 
require that the set of zero-entanglement operators has the fol- 
lowing properties ll26ll : 

1 . the set A^ is a compact and convex, such that per + (1 — 
p)(t' E Ai, whenever (T,a' € Ai and < p < 1; 

2. all matrices in A4 are Hermitian and unit-trace, so cr = 

cr^ and tr(<T) ~ 1; 

3. the set Ai includes all separable states, S C A4; 

4. the set is invariant under local unitaries, (Ua <E) 

Ub)M{Ua®UbV ^M. 




T): set of 
density 
matrices, 
cr>0 



S: separable states 



T^nAf 



FIG. 3: An heuristic diagram of different sets of zero-entanglement 
operators: operators in T> have non-negative eigenvalues; operators 
in M have partial transposes with non-negative eigenvalues; opera- 
tors in B have both the previous properties; operators in S can be 
decomposed as a positive sum of separable operators. The whole 
space shown represents all Hermitian unit-trace matrices. 



When these conditions hold, any Hermitian unit-trace matrix, 
p, can always be decomposed into zero-entanglement opera- 
tors, (0 = (1 + t)(T+ — ta~ . This result is well known for the 
set of separable states yTj, 5, and so will also hold for any 
set that contains S. Collectively, these properties ensure that 
Em is always an entanglement monotone, and similarly 12711 
for logarithmic quantities, which we denote: 



ELM{p) = ^0g2{^ + '^EM{p))- 



(10) 



Familiar monotones in this family include negativity and ro- 
bustness of entanglement, which are: 

Es{p) = Es{p), 
En{p) = Em{p), 

where S is the set of separable matrices, and Af = T)^ . The 
equivalence for robustness is immediate, but for negativity re- 
quires a little work. The more familiar definition of negativity 
is in terms of the trace norm and partial transpose. For any 
unit-trace hermitian operator, there is always decomposition 
into positive unit-trace operators p = (1 + t)(T+ — ta^ . We 
can always find such a decomposition where cr^ are orthog- 
onal, and then the trace-norm is simply 2t + 1. Clearly no 
smaller value can be achieved, and so finding the trace norm 
is equivalent to a minimization problem: 

\\p\\^mi{2t + l\p = ta+ -{l-t)(j--(j^ eV] (11) 

However, when calculating the negativity we take the partial 
transpose first: 



11/11 



inf{2i 
inf{2t 



l\p' 
1|P 



= (1 



t)a+ - ta 
-{l-t)a 



-. ± 



a' e V}, 



(12) 



Hence, when calculating the negativity, instead of taking the 
partial transpose we can find the minimal decomposition w.r.t 
the partial transpose of the physical density operators. 

In addition to more familiar entanglement monotones, Vi- 
dal and Werner also discussed an intermediate monotone for 
the set B that contains a that are positive and have a positive 
partial transpose. They observed that these monotones are re- 
lated by the inequalities: 



Es{p) > Es{p) > Em{p) > 0, 



(13) 



which follows from the inclusions S C B C Af. For an 
overview of how these sets relate see figure. |3] 



in. THE DUALITY THEOREM 

Having introduced the family of decomposition-based en- 
tanglement monotones, we now formulate our first upper- 
bound on the corresponding EC monotones. 

Theorem 1 For any logarithmic decomposition-based entan- 
glement monotone, Elm{p)> the corresponding entangling 
capacity of a CPTP map A satisfies: 

e1j^{K) < EiUiA) = snp{ELM{H^))\<y e M}. (14) 

Furthermore if AA C D, where D is the set of physical density 
matrices, then we have equality, and hence: 



Ef^iA) ^ Ef:liA). 



(15) 



This theorem instructs us that the resource assisted EC can 
never exceed the entanglement produced from applying the 
CPTP map to an operator with zero-entanglement w.rt. the 
entanglement monotone used. 

Consequently, a rough interpretation of this theorem is that 
the EC monotones considered are inherently resource inde- 
pendent. However, this is only true when the above equal- 
ity holds, as it does for log-robustness EC. However, here 
only the inequality is proven to apply for log-negativity EC. 
This is because the set of zero-entanglement operators, M, 
includes unphysical operators with negative eigenvalues (for- 
mally M % T)). Clearly, when maximization over some un- 
physical operators is required to rigorously derive an upper- 
bound, we do not guarantee that the upperbound is attainable. 

Since the resource independence interpretation is not al- 
ways strictly accurate, the more cautious reader may prefer to 
think of the theorem as simply relating two optimization prob- 
lems. The optimization giving the upperbound we will call, in 
absence of more appropriate terminology, the dual problem. 

To prove the theorem, we begin by recalling our definition 
of resource assisted entangling capacity: 

-Bf^(A) = suj>{Elm{Kp)) - Elm{p)\p € 'D}- (16) 

We shall use p* to label a physical state that achieves 
the supremum. The initial entanglement of this state is 



Elm{p*) = log2(l 
decomposition: 



2i), where there will exist an optimal 



p* = (l + <)a^ 



ta ;cr e M. 



(17) 



Think of this as a decomposition into zero-entanglement op- 
erators. Acting with the CPTP map gives: 



A{p*) = {l+t)A{a+)-tA{a-), 



(18) 



Even though A(cr^) are not always physical density matrices, 
we can still apply the formula for calculating their entangle- 
ment and we can compare it with the entanglement produced 
by acting on the physical state p* . Performing this calcula- 
tion, and some algebraic manipulation (see Appendix|C]l, we 
find: 

Elm{Kp*)) - Elm{p*) < max{£;LA4(A(a±))}. (19) 

This tells us that if we knew the optimal physical input, p*, 
then the entanglement produced would not exceed the entan- 
glement produced from the zero-entanglement operators, (t+ 
and a^. However, we don't yet know what p* is! Despite 
this, we can still give an upperbound by maximizing over all 
zero-entanglement operators, a E Ai, such that: 



^L>,(A) < snp{ELMiM'^))\<y e M}, 



(20) 



which proves the first part of our theorem. If all zero- 
entanglement operators are also physical density matrices, 
then the maximum is physically attainable, and the inequal- 
ity is saturated. Formally, if Al C P we have an equality, 
which is an especially strong result for log-robustness: 

Corollary 1 The log-robustness entangling capacity of any 
CPTP map A satisfies: 



EfsW = Els{A) 



(21) 



and so the maximum possible increase is achievable with an 
initially separable state. 

For log-negativity the result is weaker It tells us that the 
resource-assisted EC is no greater than the dual optimization 
of the maximum increase when the initial matrix is a, poten- 
tially unphysical, zero-entanglement operator, a e M. The 
sets T) and M are two very closely related species of matrices 
( note that M ~ T)^ ) and so it isn't immediately clear that we 
have gained much by exchanging a maximization over a E T> 
for a maximization over a G Af. However, the duality theo- 
rem will prove its worth by playing a pivotal role in deriving 
the entangling capacity for two-qubit unitaries. 

From a numerical perspective, the upper bound is actually 
easier to work with for reasons we explain later. For now we 
simply remark that numerical simulations indicate that little 
is lost from the relaxations needed to derive this bound, as 
numerical studies have not revealed any instances where the 
bound cannot be saturated. We discuss this further is sec- 

tionlvnl 



IV. AN UPPER BOUND ON LOG-NEGATIVITY 
ENTANGLING CAPACITY 

So far we have derived a general upper bound for a whole 
family of EC monotones. Herein we focus solely on log- 
negativity EC. Also, whereas the previous bound involved a 
maximization procedure that typically must be performed nu- 
merically, here we will derive a bound that can be evaluated by 
standard algebraic techniques. The bound we find here takes 
a very simple form for two-qubit unitaries, with later sections 
demonstrating that the bound is saturated by preparing the ap- 
propriate Choi matrix. 

We begin by defining some notation. For a CPTP map A 
there will always exist a Kraus decomposition; 



A(p) = ^if,pif|, 



(22) 



and each Kraus operator has a Schmidt decomposition: 



K 



= X! ^JV^^i.i ' 



'B. 



J, 4- 



(23) 



Recall that operator Schmidt decompositions 12811 always 
have (/) d = min(da, db)', («) coefficients A that are positive 
real numbers; and (///) sets of operators {Ai^, ..A^2 i} and 
{Bi^i, ...i?2<i.i} that are orthonormal sets, with orthonormal- 
ity defined w.rt the inner product (M, N) = tr(A/tiV). Our 
theorem will also make use of the operator norm: 



IMI 



sup 






(24) 



which for Hermitian M is simply the largest eigenvalue. We 
can now state the main result of this section: 

Theorem 2 The log-negativity entangling capacity of any 
CPTP map A is upper bounded by: 

E^^{K) < log2 f ^ ||OA.||op.||Os,«||op j , (25) 

where: 

OA,i = 2^ '^j.i^j,i^j,i 
3 

Oba = zZ^i^-'^i.i^i^'- 

3 

For unitary maps, the CPTP map has a simple form with a 
single Kraus operator The theorem simplifies even further 
for a particular class of unitaries that we call basic unitaries, 
which we define as: 

Definition 1 A unitary U acting on C"*" (E) C'^'' is said to be 
basic if there exists an operator Schmidt decomposition: 



U = J2^kAk®Bk, 



(26) 



where A^ and Bk are proportional to unitaries, such that 
AlAk = 1/da andBlBk = 1/4- 



For basic unitaries the related operator norms are simply 

I|Oa||op = ||Ob||op 

tells us that: 



J2i K/Vdadb, and so our theorem 



Ef^{U)<2logJY, 



A, 



Vdadb 



(27) 



Not all unitaries are basic, but many are. Crucially, all two- 
qubit unitaries are basic (see Appendix. IdTi. and so this result 
applies to many physical systems. 

Having laid out the results of this section, we now turn 
to proving them by utilizing the duality theorem, Thm. ([T]). 
We will find an upper bound on the resource-assisted EC, by 
considering the dual problem of maximizing the entanglement 
producible from zero-entanglement operators, a e J\f. Since 
the log-negativity is a monotonic function of the negativity, 
the maximum is achieved for the same a, so we are interested 
in finding the maximum of: 



WH'yfW = 



z^ 



iK,aK,f\l 



< 






mu<^K,f\i 



(28) 
(29) 



with the second line following by convexity of the trace norm. 
Our proof will proceed by focusing on the maximum for in- 
dividual Kraus operators. For brevity, herein we drop the i 
subscript, s.t. Ki -^ K. For these single terms we have that: 



(30) 



A'crA't = Y^ XkXjiAk ® Bk)a{Aj ® Bj)\ 

k,3 

taking the partial transpose we arrive at: 

{KaK^f = Y,^kXMk®Bj)cj^{Aj®Bl)\(3\) 

k,j 

The next step is to find the trace norm of this expression. As 
outlined in section |II] taking the trace norm is equivalent to 
a minimization problem over all decompositions into positive 
unit-trace operators. By finding a decomposition, which does 
not necessarily give the minimum, we are able to deduce an 
upper bound on the trace norm. Here, we are always able 
to find a decomposition (see Appendix |E1i such that the trace 
norm is bounded as follows: 

WiKaK^fW < Y^XkXMA^^i ® BlBk<y^)- (32) 
j.k 

Making use of the more compact notation introduced earlier, 
this is equivalent to: 

Maximizing over all cr e A/" (or cr'" e 2? ) is equivalent to 
finding the operator norm, and so: 

sup||(i^aA't)r|| < \\0a®01\\o^, (33) 



aeN 



fc=i 



||Oa||op-||ObI|op, 
||Oa||op-||0_b||op. 



The Hilbert space on which 
a Hermitian map A acts. 




An ancillary Hilbert 
space of equal dimension 



FIG. 4: The Hilbert spaces used in defining tlie Clioi matrix pA ~ 
(A ® Il)($) for a dynamical map, A. For any bipartition, Ha ® Hi,, 
of the space on which A acts, we find the Choi matrix can be defined 
with the initial state as a separable state, |$) = |$a)|$s), w.r.t this 
bipartition. 



The second line follows because the operator norm is mul- 
tiplicative under tensor products, and the last line follows be- 
cause transposition does not change the eigenvalues of Hermi- 
tian matrices. Taking the logarithm and applying the duality 
theorem (Thm.[T]), proves theorem|2] 

It is important that this result holds when local ancilla are 
available, and we will clarify this now. For brevity, we have 
denoted the Kraus operator by K. However, we implied the 
availability of local ancilla, which is equivalent to finding the 
EC for an operator i^ (X) Ha'^^b'- For this extended Kraus op- 
erator, the appropriate bound is determined by \\0a <E) la' I lop 
and 1 1 Ob® lfc'||op. However, ||Oa (81 la'||op = ||Oa||op, and 
so the upper-bound does indeed allow for local ancilla. 

As we remarked earlier, this result takes an especially sim- 
ple form for basic unitaries, and we shall see that basic uni- 
taries saturate the bound. However, in section IVIIII we give 
some explicit examples that fail to saturate it. 



V. CHOI MATRICES 

Before showing that the bound can be saturated, we will re- 
view the concept of a Choi matrix. The Choi-Jamiolkowski 
isomorphism jig, llTT shows that there exists a one-to-one 
correspondence, an isomorphism, between Hermitian matri- 
ces and Hermitian maps ll29ll . In particular, for every map A 
that acts on a Hilbert space, C^" (^ C''^ there exists a Hermi- 
tian matrix p\ that acts on C^" (E) C'''' . We call pA the Choi 
matrix for A, and it is defined by: 



PA == (A® ]!)($), 
where <1> = |$)($|, and: 



dadb 



1$) = ^ \x)a.b\x)a',b'/\/dadb 



(34) 



(35) 



This state is entangled between Hilbert spaces Ha^b and 
Ha'.b'- However, the initial state, $, can be separable with 
respective to a partition between Ha. a' and i/h^bs by choosing 
the basis such that \x)a,b = \j)a\k)b and \x)a'.b' = \j)a'\k)b' 
(see figurelH). It then follows that: 

da dh 



\<^) =Y.Y.\^'^-\^)b\j)a'\k)b' / Vd^b- (36) 

J=l k=l 

Switching the order of the qubits from a, &, a', h' to a, a', 6, h' , 
we arrive at: 



da db 



I*) = T.J2\3)-\-})-'\'')b\k)b'/VdJ'b, 

j=l fe=l 



(37) 




which is simply a separable state: 



(38) 



As such the Choi matrix is both mathematically intriguing and 
operationally meaningful as it can be prepared from a separa- 
ble state and a single use of a dynamical map. 



VI. CHOI MATRIX ENTANGLEMENT 

Now we show that basic unitaries saturate the upper-bound 
described by theorem (|2]), and more concisely by equa- 
tion dZTJ ). We achieve this by using a basic unitary to repro- 
duce its Choi matrix. Operationally, the Choi matrix is pro- 
duced by applying U to the initial separable state, |$a) I'&s)- 
Establishing the entanglement of the Choi matrix doesn't ac- 
tually rely on the operator being a basic unitary, and can be 
considered as a useful lower bound for the EC of any general 
Kraus operator Specifically, we show that: 

Theorem 3 The log-negativity entangling capacity of any 
Kraus operator K satisfies: 



Ef^{K) > Eln{K) > 21og2 ( ^ -^ 



(39) 



where Xj are the Schmidt coefficients of K. Equality holds if 
K is a basic unitary, such as a two-qubit unitary. 

It is important to remember that we can produce the Choi state 
from an initially separable state, and hence are proving that the 
log-negativity is resource independent for all basic unitaries. 
We begin with: 



K\^a)\^b) 



k 



Xk{Ak\<^A))®{Bk\^B)), 



(40) 



x=l 



and we shall prove our result by showing that, except for a 
factor of ^/dadi, this decomposition is aheady a Schmidt de- 
composition. We proceed by proving a straightforward, and 
possibly well-known, fact: 



Lemma 1 For orthonormal sets of operators {Ak\, the set of 
states llV'Afc) = Vd^Ak\^ a)} ore orthonormal. 

The lemma is proven in Appendix IF] From this lemma, it 
follows that a Schmidt decomposition of the Choi matrix is: 

Afe 



b 



K\^a)\^b)=Y. 



\/dadb 



I^AkMi 



(41) 



It is well known that the log-negativity of a pure state follows 
directly from its Schmidt coefficients, such that: 



E 



LN 



(i^$A-t) =. 21og2^ 



A, 



Vdadb 



(42) 



which holds for any Kraus operator. For basic unitaries, this 
achievable entanglement and the upper bound coincide, giv- 
ing a closed form for the EC, and proving our main result 
(theorem[3]). 

We have shown that for a particular class of unitary oper- 
ations, which includes all two-qubit gates, the maximum in- 
crease of log-negativity from one application of the unitary 
has a simple closed form. Furthermore, the same choice of ini- 
tial state produces this maximum. As such, if we are promised 
that a unitary belongs to the class of basic unitaries, then no 
further analysis is required to find the optimal strategy for gen- 
erating entanglement. We also have the surprising result that 
unlike entropic quantities, such as von Neumann entropy, the 
log-negativty EC is resource independent for 2-qubit unitaries. 

Although the Choi matrix entanglement equals the log- 
negativity EC, this does not entail the same for the EC based 
on different monotones. Indeed, Kraus and Cirac (sll showed 
that preparing the Choi matrix does not always maximize 
the linear-entropy (assuming pure states only). The counter- 
example they gave employed unitaries of the form U = 
exp'^^, where W is the two-qubit swap gate, and the Choi 
matrix was non-optimal when 9 was below a critical value. In 
these cases, there is no unique answer to what the best strategy 
is, as both entanglement monotones are meaningful. However, 
the strength of our result is that it provides a nice analytic so- 
lution, whereas no general solution for the linear-entropy EC 
is yet known. 



VII. NUMERICAL EVIDENCE 

What do the results so far tell us about the existence of re- 
source independent EC monotones? To recap, we know that 
log-robustness EC is always resource independent, and that 
log-negativity EC is resource independent for basic unitaries, 
such as two-qubit gates. The most obvious open question is 
whether log-negativity EC is resource independent for all uni- 
taries, or indeed all CPTP maps, or whether this is a singu- 
lar feature of 2-qubit gates. Here we report numerical evi- 
dence in favor of the conjecture that for all unitaries we have 
E]^^iU) = Eln(.U). 

We consider 50 different 2-qutrit unitaries randomly se- 
lected from the Haar measure, which are typically not ba- 
sic unitaries. To calculate Elis!{U) we performed a Monte- 
Carlo search over the set of pure separable states. We can 
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FIG. 5: Numerically calculated log-negativity EC for randomly se- 
lected two-qutrit gates, allowing for two-qutrit ancillary qubits. The 
achievable log-negativity with separable states (resource unassisted), 
Eln{U), is compared to EY^f^{U), which is an upper bound on 
the maximum resource unassisted increase E'£j^{U). For all exam- 
ined instances these quantities are equal (up-to numerical precision 
of lO"*^). 



ignore mixed states because the maximum for any convex 
function over a convex set is attained by an extremal mem- 
ber of the convex set. We calculated an upper bound of the 
resource-assisted EC by considering the dual problem of max- 
imizing over zero-entanglement operators. Again this opti- 
mization can be performed over just the extremal members 
of the zero-entanglement set M. Alternatively we could have 
directly calculated the resource assisted entangling capacity, 
E]^j^{U), but this would have required llsoll searching over a 
larger space including non-extremal points. Surprisingly, this 
more comprehensive search was not necessary as the upper- 
bound proved to be tight in all observed instances, as can be 
seen in figure |5] 

We allowed for ancilla qubits equal in dimension to the tar- 
get qubits, and so considered 4-qutrits in total. Since local 
maxima are not guaranteed to be global maxima, and the di- 
mensionality {3'^ = 81) involved in the problem is quite large, 
a fairly intensive search must be performed to encourage con- 
fidence in having found the global maxima. Numerical ac- 
curacy was benchmarked against unitaries randomly selected 
from a family of basic unitaries, for which the results always 
matched our analytic result up-to our working numerical pre- 
cision of 10^®. Similarly for non-basic unitaries, all instances 



of numerically found values for E^j^{U) and Elis[{U) dif- 
fered by no more than this level of precision. 

Consequently, it is extremely likely that all two-qutrit uni- 
taries are resource independent. In other aspects of entangle- 
ment theory where special properties have held for two-qubits, 
but not for higher dimensions, the emergent behavior can be 
observed for just two-qutrits. For example, bound entangle- 
ment 1I20I1 an d in-equivalence of entangling and disentangling 
capacity 12311 both manifest themselves for qutrit systems and 
larger. As such there is little reason to suspect that two-qutrit 
gates are exceptional in any way, and we are inclined to con- 
jecture that resource independence wiU hold true for all two- 
qudit unitaries. 



Vin. TRACTABLE QUDIT EXAMPLES 



IX. CONCLUSIONS 



We have seen that various analytic bounds on EC are satu- 
rated for two-qubit unitaries, and that the duality upper bound 
-^L7v(^) seems to be tight in higher dimensions. In this sec- 
tion we examine some analytically tractable two-qudit uni- 
taries, and show that in higher dimensions (/) preparing the 
Choi matrix is not always optimal; and (//) the upper-bound 
expressed by theorem |2] is not always saturated. There are 
many ways to generalize two-qubit control-not gates, with 
some generalizations giving basic unitaries, but we shall con- 
sider the non-basic (for d > 2) family that acts on systems a 
and b of dimension d: 



(43) 



C/d = (l-|0)(0|)®]l + |0)(0|®Xd, 
where Xj, is the generahzed Pauli bit-flip operator: 

d-l 



(44) 



with ® indicating addition modulo d. A slight rearrangement 
gives the Schmidt decomposition with correctly normalized 
operators: 



Ud = Vdid-l) 



1 



V^R^l) 



7^^^'«^<«'^^- 



Given that the Schmidt coefficients are \Jd[d — 1) and ^/d, 
the entanglement of the Choi matrix is simply (see Eqn. 



ELNiUd'^Ul) - 21og2 ({yjd{d-l) + ^d)/d) , (45) 

which vanishes in the limit of large d. 

However, if we consider the input state jf/') = (|0) + 
|l))|0)/-\/2, then it is easy to check that C/d|V') produces one 
e-bit of entanglement, or log2(2) of log-negativity. Conse- 
quently, for d > 2, preparing the Choi matrix is not optimal. 
No operator with only two non-zero Schmidt coefficients can 
produce more than log2(2) log-negativity and so this must be 
the entangling capacity for all values of d. 

We now consider the value of the upper-bound imposed by 
theorem|2] which tells us that: 

ELN{Ud) < log2(||OA||op.||Oi3||op), (46) 

= l0g2 ((Vd) . (Vl-rf-1+Vd^)) , 

= logsfVd^T+i), 



which goes to infinity as d increases, and hence this is also not 
tight against the actual EC of log2(2). 

On a more positive note, we observe that these unitaries are 
also resource independent. Again this provides further support 
for our conjecture. 



We have investigated the maximum entanglement that can 
be achieved by a dynamic process allowing for use of ancil- 
lary Hilbert space, with particular attention paid to the boost 
that can sometimes be achieved by exploiting prior entangle- 
ment. These quantities depend on the underlying entangle- 
ment monotone that is used to quantify the entanglement of 
the quantum state. We considered a family of monotones 
that includes log-negativity and log-robustness, with a par- 
ticular focus on log-negativity. We found that this family 
of entangling capacity monotones is often insensitive to the 
availability of prior entanglement, a property we call resource 
independence of the entangling capacity. Our most general 
result, which covers a whole family of entangling capacity 
monotones, we call the duality theorem as it tells us that the 
resource-assisted entangling capacity cannot exceed the so- 
lution of a dual problem determined by the structure of the 
entanglement monotone. 

The duality theorem immediately entails that log- 
robustness is resource independent for all dynamical opera- 
tions. For log-negativity, the duality theorem entails a weaker 
result with a murkier interpretation. However, further results 
built on the duality theorem prove that log-negativity entan- 
gling capacity is resource independent for all two-qubit uni- 
taries, and other so-called basic unitaries. Furthermore, we 
presented numerical evidence that two-qutrit unitaries are re- 
source independent. This leads us to the conjecture that all 
unitary gates are resource independent for this particular met- 
ric of entangling capacity. Settling this conjecture is an espe- 
cially interesting open problem. 

We have also shown that the log-negativity entangling ca- 
pacity of two-qubit gates has a closed form, and equals the 
entanglement of the corresponding Choi matrix. For no other 
measure of entangling capacity has a closed form been found, 
and so this result may prove extremely useful in simplifying 
the analysis of entangling capacity. We have seen that beyond 
two-qubit unitaries the Choi matrix no longer captures the log- 
negativity entangling capacity. However, the simplicity of the 
result for two-qubit unitaries is encouraging evidence that a 
closed form may be also be found for higher dimensional uni- 
taries, and maybe even general dynamical maps. 
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Appendix A 



Appendix B 



In the papers preamble we asserted that subadditivity of an 
entanglement monotone is a necessary property for the corre- 
sponding EC to be resource independent. The basic intuition 
is that when subadditivity is violated, one can boost entan- 
gling capacity by simply appending an ancillary entangled 
state that actually plays no active role. Indeed, this spuri- 
ous boost may ensure an infinite amount of resource assisted 
EC for almost all CPTP maps, making the quantity devoid of 
concrete meaning. This rough intuition will suffice for un- 
derstanding our main results, but for the interested reader we 
shall provide a rigorous proof here. 

If an entanglement monotone is not subadditive, then there 
must exist at least two quantum states pi and p2, such that: 



E{pi®P2) >E{pi) + E{p2). 



(Al) 



We now define a CPTP map that will violate resource inde- 
pendence: 

Ai(a) = tra,fc(|0,0)(0,0|a)pi + (l-tr,,b(|0,0)(0,0|a))l/d, 

(A2) 
where d is the dimension of pi, and trj,, (,(...) is a partial trace 
over target qubits a and h. Clearly, with the initially separa- 
ble state |0, 0) and one application of Ai we can produce the 
output pi. For any other separable input, ga® crs, we have 
that: 

Ki{aA®cfB) = ?7Pi(8)cra'(g)o-b' + (l-?7)-®cr^,«)cr[,,, (A3) 

a 

where: 



We show here that an entanglement monotone, Ei, that 
is subadditive, can always be rescaled into a non-subadditive 
monotone. Consequently, a resource independent EC can al- 
ways lose this property by appropriate rescaling. A second 
monotone £'2 is a rescaling of Ei if and only if there exists 
a monotonically increasing function, g, such that E2{p) = 
g{Ei{p)). We proceed under the weak assumption that if 
Ei{p) > 0, then 



Ei{p®p) > Ei{p). 



(Bl) 



This assumption is satisfied for all well-known entanglement 
monotones, and only fails to hold for extremely contrived con- 
structions. We can now break subadditivity for p (g) p, by 
choosing a rescaling function such that: 



gix + S)> 2g{x), 



(B2) 



where x = Ei{p), and x + 5 = Ei{p ® p). Our earlier 
assumption ensures that 6 > 0, and so we simply choose a 
rescaling function that more than doubles in the gap between 
X and X + S. 

Since a rescaling can always break resource independence, 
it is reasonable to ask whether the inverse operation can al- 
ways be performed. That is, can we always rescale a mono- 
tone to create a resource independent EC? We do not cur- 
rently have a solution to this problem. However, one fact is 
clear: rescaling a strongly-additive resource-dependent entan- 
glement monotone can not make it resource independent with- 
out also breaking strong additivity. 



77 = tr(|0,0)(0,0|a), 

aa'(E)<Jb' = tra,b{\0,0){0,0\(7)/T], 

aa-®<Jb' = tra,fc((l-|0,0)(0,0|)a)/(l-77). 

The above output state can always be produced from pi by lo- 
cal operations and classical communication. Therefore, since 
E is monotonic under such operations the resource unassisted 
EC cannot exceed E{pi), and we have E{Ai) = E{pi). 

We will now prove that the resource assisted EC can al- 
ways exceed this value. We consider the entangled input state 
|0,0)(0, 0| (g) p2 where p2 is in the ancillary Hilbert space 
Ha' (E> Hf. Since the ancillary qubits are uncorrected w.rt 
the target qubits, this ancillary resource plays no active role in 
the dynamical process, such that: 



Ai(|0,0)(0,0|)®p2=pi®p2, 



(A4) 



the entanglement of which is: 



£;(Ai(|0,0)(0,0|)®p2) = E{pi<E>p2), (A5) 

> E{p,) + E{p2), 

where the second line uses the violation of subadditivity. It 
follows direcdy that i?'^(Ai) > E{Ai), and so the EC is not 
resource independence. 



Appendix C 

In this appendix we compare the entanglement of A(p*) 
with that of A(o'^), which results in the inequality in Eqn.[T9l 
Recall that the zero-entanglement operators, a^, are given by 
the decomposition of p* in Eqn.fTTl 

We begin by considering A{a^), where the CPTP map, 
A, preserves the trace and Hermiticity of a^, and hence the 
new operators also have optimal decompositions into zero- 
entanglement operators, such that: 



A(cr^ 

A{a- 



(1 
(1 



- a)a 
b)a- 



aa 
ba- 



(Cl) 
(C2) 



Combining these equations with Eqn. [18] we have that: 



A(p*) = {l+t)il + a)cr+'+ -{l+t)aa+'-, (C3) 
-t{l + b)a^'+ +tba^'^. 



We can simpHfy this expression by defining new operators 

{I + t){l + a)a+-+ + tba--- 



a' = 



a^ 



{I + t){l + a) + tb 
{l + t)aa+^- +t{l + b)(7-'+ 
(l+t)a + f(l + 6) ■ 



(C4) 
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Since the zero-entanglement operators form a convex set, it 
follows that CT^'^ are also zero-entanglement operators. In 
terms of these new operators we have: 



l\~l 4-1-2 



where. 



A(p*) = {l + t'Ya^ -t'a 



i' = (l+i)a + i(l + 6). 



(C5) 



This decomposition of A(/9*) is of the correct form for cal- 
culating the entanglement, but may not be the decomposition 
that minimizes t' . Hence, it gives an upper bound on the en- 
tanglement: 

Em{Hp*)) < t' = {{l + t)a + t{l + b)). 

Next, we define c = max{a, 6}, and use it to replace these 
variables. This replacement further increases the larger side 
of the inequality, and so we have that: 

^>t(A(p*)) < (l + t)c + t(l + c). (C6) 

Adding contributions to both sides gives a useful factorization 

1 + 2Em{Kp*)) < (l + 2c)(l + 2t), (C7) 

= (l + 2c)(l + 2^^(p*)). 

Taking the logarithm, and applying Eqn (fTOl i, gives: 

Elm{Hp*)) < log2(l + 2c) + Elm{p*)- (C8) 

Subtracting the initial entanglement gives: 

Elm{Mp*)) - Elm{p*) < log2(l + 2c). (C9) 

Notice how the terms involving the variable t have cancelled 
out, simplifying the nature of the bound. Next, we observe 
that the following three facts (/) c = maxja, 6}; (//) 1 + 2a = 
^LA4(A(cr+)); and(/;7) 1 + 26= £'LA4(A(cr~));joinly entail 
that: 



l + 2c = max{^L^(A(CT±))}, 



(CIO) 



and hence: 

Elm{Mp*)) - Elm{p*) < max{£;LA4(A(a±))}. (Cll) 
This complete derivation of the inequality. 

Appendix D 

In general, to search for a basic decomposition of a given 
U we find an orthonormal Schmidt decomposition. This will 
be a basic decomposition if each term is proportional to a 
unitary. If all the non-zero Schmidt coefficients are differ- 
ent, then there is a unique Schmidt decomposition of U, and 
finding it gives a conclusive answer as to whether U is basic. 
Whereas degeneracy in the Schmidt coefficients entails a fam- 
ily of valid Schmidt decompositions and the problem is more 



involved; in these instances one must exhaustively search the 
family of Schmidt decompositions for one that is basic. 

However two-qubit unitaries are always basic and we shall 
give a proof here. This feature of two-qubit unitaries follows 
from the widely exploited I4i, [8i, i9i] fact that all two-qubit uni- 
taires are local unitary equivalent to a diagonal operator f/diag, 
such that: 

U={UA®UB)Udi.g{VA®VB), (Dl) 

where the diagonal operator has the form: 

C/diag = exp'^^U-^-^^-^)^ (D2) 

where 0-1.2,3 are the Pauli spin matrices and the operator is 
unitary for all real Cj . Expanding out the matrix exponential, 
and including the local unitaries we have: 

4 
U={Ua<» Ub){Y^ fljCTj ® crj){VA ® Vb), (D3) 
3=0 

where gj are complex numbers determined by Cj . We define 
Aj = UACTjVA/Vd^ and B-j = Aig{aj)UBcrjVB/Vdb, with 
any complex phase absorbed by the operators. It is straight- 
forward to check that {Aj} and {Bj} are orthonormal sets 
and are proportional to unitaries. Consequently, all two-qubit 
unitaries have a basic decomposition with Xj = \aj\. 



Appendix E 

This appendix takes equation ( |3TI ). and derives an upper 
bound on its trace norm, producing equation ( l32b . As outlined 
in the main text, the overall strategy is to find a valid, but not 
necessarily minimal, decomposition. We begin by introducing 
a more compact notation: 



(KaK^) 



r\r 






where. 



Afj, fe = {Ak ® Bj)<7''{Aj ® Bj)t. 



(El) 



(E2) 



The trace norm is evaluted by searching over real decomposi- 
tions of positive unit-trace operators. Recall that for any posi- 
tive Hermitian matrix, p, the matrix CpC'' is also positive and 
Hermitian for any C. Consequently, the diagonal terms Mk,k 
are always positive and Hermitian. However, the cross terms 
Mj^k^j are not individually positive, and so the trick is to find 
an appropriate decomposition for collections of cross terms. 
It turns out that we can find such a decomposition for the sum 
of paired cross terms, Mj^k + M^.j, such that: 



M. 



j,k 



Ml 



k,J 



M. 



j,k 



^' 



(E3) 



where the new operators are: 

Mf:^ = l{X,,k±Xk,,V{X,,k±Xk,-)\ (E4) 
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and the Xj^k are operators: 



X 



3,k 



A, 



'Bj. 



(E5) 



Now our tilded operators, il/n., are positive matrices. Col- 
lecting these equations together, we have that: 

k k.j<k 

where each of the operators is positive and Hermitian. Now 
we can use our decomposition, which is not necessarily opti- 
mal, to prove a bound on the negativity, such that: 

k 
k,j<k 

Again we can further simplify pairs of terms: 



tr(M+,+A/rJ 



= tT{Xj^k^^xl^+Xk,ja^xl^), (E6) 
= tT{xl^X^^k^^)+tT{xl^Xk,ja^), 



where the second line follows from the linearity and cycli- 
cality of the trace. Since the diagonal terms similarly satisfy 
tr(Mfcfc) = tT{Xl i^Xk^a^), we have that: 

3,k 
■J,k 



This completes our proof. 



Appendix F 



Here we prove lemma. [T] as follows: 



{^^aMa,) = da{^A\{Al®l){A,®l)\^A) (Fl) 

= ^(.T,a:|(4®l)(A,®l)|y,2/) 



x,V 



Y,{x\AlA,\y){x\y) 
Y.{x\AlA,\x). 



The final line is simply the trace operation performed in the 
computational basis, and so: 



(V'aJV'a,) = tr(4A,) 



(F2) 



3fcj 



where we have used that the operators {Aj} are orthonormal. 
A similar proof holds for system B. 
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